General Theorem about Gravitational Lensing 
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ABSTRACT 



We study the general theorem about gravitational lensing which states the 
relationship between the numbers of images with different parities. Our formulation 



allows an extension to the nontransparent and singular model. 
Subject headings: gravitational lensing 
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There are the general theorems about the gravitational lensing( MacKenzie 1985| ; Schneider, 



Ehlcrs, &: Falco 1992J ) which mainly state the relationships between the numbers of the images 



with different parities; 

V '. n i ~ n u + niu = 1 , (1) 

where nj is the number of images of type I. We use njj and nm analogously. 

The types of images are defined by the Jacobian matrix of the Fermat's potential (f>, 

-^aSfc (! ' J = 1 ' 2) ' (2) 

Type I : both eigenvalues of the Jacobian matrix are positive. 
Type II : one eigenvalue is positive and another negative. 
Type III : both are negative. 
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Therefore type I and type III images have positive parities and type II images have negative ones. 

The existing general theorems are general for transparent and non-singular lens models. 
Typical examples are the elliptical lens models with softened power law behavior (Binney & 



|Tremaine 1987 ; Schneider, et al. 1992). However, they do not exhaust all the lens models. 



They are rather phenomenological models, whereas there is a widely valid model. It is multipole 
expansion model(MPE model) which is obtained by expanding deflection potential with respect 
to impact parameter ( [Yoshida fc Omote 1988 ; |Kakigi, Okamura, &: Fukuyama 1995 ; Fukuyama, 



Kakigi, &: Okamura 1997| ). The latter is less popular than the former in lensing analysis but is 



very powerful (Fukuyama, et al. 1997). The MPE model is nontransparent and singular, and we 
see how the general theorems are extended to incorporate MPE model and other models. In this 
extension we propose the generalized index theorem valid also to the nontransparent and singular 
lens model and the formulation in complex coordinates is found to be very convenient. 



2. Formulation 

Before proceeding to the concrete physical situation, let us begin with the mathematical 
preparations. Index theorem for critical points in two dimensional flat or closed curved surfaces 
is well known for mathematicians ( Eguchi, Gilkey, fe Hanson 198C| ). They are summarized as 



follows: We consider the vector field on the surface S. Denoting the critical points (at which the 
vector field vanishes) by Pi on S, we get the following theorem. 

/^ (Index of the vector at Pi) = x(S) , (3) 

i 

where x(S) 1S the Euler characteristic of S. For flat space x is one an d the result becomes equation 
(Jl]). This theorem does not hold in the case that the vector field has singular points. So we 
reformulate the index theorem so as to incorporate pole terms. 

Let us consider a complex function f(z) with a complex variable z that is not necessarily 
holomorphic. The function / plays the role of the vector field mentioned above. The index of / is 
defined by 

Ind(/) = ^r£ din (/(*))■ (4) 

This measures the variation of the phase of / along the closed curve C z which goes anti-clockwise. 
When C z does not enclose any zero points or poles, the index vanishes. If C z encloses some zero 
points or poles, {Pi}, then the index is given by the summation of the contributions from each 
point Pi, 

Ind(/) = £lnd(/;P 4 ), (5) 

i 

Ind( /; P,) = ±-J c(i) d In (f(z)), (6) 



(i) 
where C z is the circle enclosing the point Pi. In the context of the gravitational lenses, / is 

replaced by 4>z(= dcj>/dz; z is complex conjugate of z) and 4> z = is just lens equation (see eqs.[14| 

and |15[|). Therefore ordinary images are zero points of order one. If lens object has finite total 

mass and extended regular mass density, then we have the pole only at infinity. Whereas in the 

case of MPE models and models with singular mass density such as point masses, other poles are 

produced. 

From the form of equation (||) we should pay attention to each contribution. Let us consider 
the map u = f(z) in the neighborhood of Pi. Then the Ind(/ ; Pi) is rewritten as 

Ind(f ;Pi) = ^-f dlnu, (7) 

where C w = u){C z )• The value of the right hand side is the (algebraic) number of times that we 
wind around C^ when we go around C z in the anti-clockwise direction. 

Firstly we consider zero point case and proceed to the case of poles. Hereafter we set the 
origin of complex z-plane at the zero point Pi and we use bold-face letters for the value of any 
function g, for example, g = g(Pi). 



2.1. The case of zero point of order one 

In the neighborhood of a zero point, uj can be written as 

u = zi z + zi- z + 0(\z\ 2 ) , (8) 

where lower indices z and z mean partial derivatives with respect to them, for example, 

Qk+l 

fzk2 ' ~ dz*d¥ f ' 

We define fi by 

h 

Then the Jacobian D of the mapping from z to w is in this case given by 

Z?=|f z | 2 -|f 5 | 2 = |f 2 | 2 (l-| / u| 2 )- (9) 

That is, u) = f(z) reserve (reverse) the direction of rotation if |//| < 1 ( |/i| > 1 ). 
Winding number of C u is one for the point of order one and 

M/ifl)-/ 1 , ! or j"j^ do) 

! —1 tor \fj,\ > 1 . 



2.2. The case of zero point of higher orders 

About a zero point, u is expanded as a polynomial of z and z. However the winding number 
is determined by the lowest power term as is easily checked. 

For the simplicity, we consider the case that dominant part in the neighborhood of the zero 
point consists of the only term as 

u> = z z z z + (subdominant terms) (11) 

k\ l\ 

We call this type of zero point as (k, l)-type critical point for later use. Although equation (|ll|) 
restricts the form of /, as we see later, for the case of point masses and MPE models, it is enough 
to examine equation (|TT|). 



Then the index of / is given by 

lnd(f ;P i ) = k-l. (12) 



Therefore, from the previous subsection |2.1| , we can regard zero points with |/x| < 1 as (l,0)-type 
and ones with \fi\ > 1 as (0, l)-type. 



2.3. The case of pole terras 

As was mentioned, besides the zero points, pole terms of / also contribute to Ind(/). In the 
neighborhood of pole of / we define g by g = 1//, then 

lnd(f;P i ) = -lnd(g;P i ) (13) 

comes from the definition of equation (0). If / consists of several pole terms, Ind(/) is determined 
by the highest pole term. 

Thus, similarly to the case of zero points, we can regard the pole whose dominant term has 
the form of equation (O) as (fc, Z)-type critical point with negative integers k and I. 



3. Application to gravitational lensing 

We have completed mathematical preparations and proceed to application to gravitational 
lensing. In this case / is defined through the Fermat's potential by 

/ = 4>z , (14) 

where 

ct> = \\z-z s \ 2 -^{z). (15) 



z(z s ) are image (source) positions in the complex notations and ip is the deflection potential 
defined by 

ijj{z) = - f d 2 z'K(z') In \z - z'\ , (16) 

vr J 

where the integral measure d 2 z means dSR.(z)d$s(z) and k is normalized surface mass density of 
lens object. We assume 



lim 

\z\— *oo 



0, (17) 



z 
and this condition may hold for the models with finite total mass. 

For the large impact parameter, (f) s becomes 

(j>z = \{z-z s ) + o{\z\) . (18) 

That is, infinite point becomes pole. Remarking that the anti-clockwise loop C z (\z\ S> 1) is 
viewed as clockwise one from the infinite point, we obtain the contribution of the infinite point to 
the index as 

Ind(cfe) = 1 . (19) 

Diminishing the radius of the circle C z , it is entangled with the zero points and poles. Then 



equation (19) must be equal to the sum of these contributions. 



As we have mentioned, the zero point of order one in the neighborhood of the origin is 
classified into (l,0)-type for |/z| < 1 and (0, l)-type for \/j,\ > 1 by the contribution to the index. 
Thus, we can say that the contribution comes from the term 

Ck,l k -i 

where C^i is a constant and positive (negative) k and I correspond to the zero points (poles). We 
denote the number of (k, Z)-type critical points as n' ' '. Then from the arguments of the section 
|, we get 

1 = Ind(^) = J2( k ~ n(M) • ( 2 °) 

k,l 

• For the extended lens models which have no pole we obtain 

1 = ni 1 * ) - n^- 1 ) = n+ - n_ = m + n HI - n H , (21) 

where n + (n^) is the number of the images with the positive(negative) parity. 

• For the point masses, we have a pole of (0, — l)-type for every point mass and therefore 

1 = m + nni — nn + (# of masses) . (22) 
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• For a MPE model up to 2 fc -pole, Fermat's potential is given by(Kakigi, et al. 1995) 

^=>-^| 2 -^n(--) + Efe + ^) (23) 

n=l \ / 

and the vector field becomes 

So there give rise a pole of (0, —k — l)-type and 

1 = n (1 ' 0) - n^°' l) + (Jfe + 1) n^- k -^ = m + n IH -n n + k + l. (25) 

Lastly, we give the concrete example of the theorem in the multiple images of 
PG1115+080 (|Young, et al. 198 1| ; [Christian, Crabtree, fe Waddell 1987J ; [Kristian, et~al 



1993 ; |5chechter, et al. 19971 ). Fig- H] shows critical and caustics lines of the Fermat's potential 



and image positions numerically fitted by a model with softened power law surface density ( Asanoj 



fc Fukuyama 1996|). Critical lines are those on which Jacobian determinant equation (H) vanishes. 



So the parities of images are reversed across a critical line. Images outside the outermost critical 
line have positive parities. Hence critical lines are useful for determination of parities of images. 

From positions of the critical lines, we can know that A\ and C belong to nj, A<i and B do to 
nn and D to nm. So the calculated numbers of images satisfy the original theorem equation (|l]). 

On the other hand Fig. ^ shows fitting of the same lens system by a MPE model including up 
to quadrupole (Fukuyama, et al. 1997). A2 and B belong to nj. The other A\, C, D and E do to 
nn. Thus, these satisfy equation (p5[) for the case of quadrupole, k = 2. 



The authors would like to thank Professor E.L. Wright for comments. 
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Fig. 1. — The observed image positions (crosses) and numerically fitted ones (diamonds) by a model 
with softened power law surface mass density. The square shows source position. The solid and 
dotted lines are critical line and caustics, respectively. 
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Fig. 2.— ibid, by a MPE model. 



